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ABSTRACT 

We investigate the formation of polarization profiles induced by a magnetic 
field in the He I multiplet at 10830 A. Our analysis considers the Zeeman splitting 
in the incomplete Paschen-Back regime. The effects turn out to be important and 
produce measurable signatures on the profiles, even for fields significantly weaker 
than the level-crossing field (~400 G). When compared to profiles calculated with 
the usual linear Zeeman effect, the incomplete Paschen-Back profiles exhibit the 
following conspicuous differences: a) a non-Doppler blueshift of the Stokes V 
zero-crossing wavelength of the blue component; b) area and peak asymmetries, 
even in the absence of velocity and magnetic gradients; c) a ~25% reduction in 
the amplitude of the red component. These features do not vanish in the weak 
field limit. The spectral signatures that we analyze in this paper may be found 
in previous observations published in the literature. 

Subject headings: line: profiles - Sun: atmosphere - Sun: magnetic fields - Sun: 
chromosphere 
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1. Introduction 

Many of the physical challenges of so- 
lar and stellar physics arise from magnetic 
processes taking place in the outer regions 
of the atmosphere. Unfortunately, our 
present empirical knowledge on the mag- 
netism of the chromosphere, transition re- 
gion and corona is still very crude. In our 
opinion, there are two key avenues that 
should be pursued to improve the situa- 
tion. 

Firstly, new instrumentation is needed, 
including a space-borne UV polarimeter 
and a large ground-based solar telescope 
optimized for (spectro-)polarimetric obser- 
vations (e.g., the Advanced Technology So- 
lar Telescope). Secondly, but equally im- 
portantly, we need appropriate diagnostic 
tools that will permit us to infer the mag- 
netic field vector from the observed polar- 
ization in suitably chosen spectral (chro- 
mosphcric and coronal) lines. The present 
paper represents a contribution to this sec- 
ond line of research. 

In particular, our aim here is to demon- 
strate that the determination of the mag- 
netic field vector via the analysis of the ob- 
served polarization in the He I 10830 A 
multiplet must be carried out considering 
the wavelength positions and the strengths 
of the Zeeman components in the incom- 
plete Paschcn-Back effect regime. Previ- 
ous research done with this multiplet has 
ignored this effect (with the only exception 
of Trujillo Bueno et al. 2002), which re- 
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ence Foundation. 
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suited in considerable confusion surround- 
ing the formation physics of these lines. 

The He I 10830 A multiplet originates 
between a lower term (2^Si) and an upper 
term (2^P2,i,o)- Therefore, it comprises 
three spectral hues (see, e.g., Radzig & 
Smirnov 1985): a 'blue' component at 
10829.09 A with J/ = 1 and J„ = 
(hereafter Transition 1, or Trl for ab- 
breviation), and two 'red' components at 
10830.25 A with J„ = 1 (hereafter, Tr2) 
and at 10830.34 A with J„ = 2 (here- 
after, Tr3) which appear blended at solar 
atmospheric temperatures. As shown in 
Fig 1, the energy difference between the 
levels with J„ = 2 and J„ = 1 is such that 
their respective magnetic substates cross 
each other for magnetic strengths between 
400 G and 1500 G, approximately. There- 
fore, at least in this range of strengths, the 
Zeeman splitting of such J-levels is compa- 
rable to their energy separation. As a re- 
sult, the perturbation theory of the famil- 
iar Zeeman effect is no longer valid. The 
energy positions of the magnetic substates 
as a function of the field strength have to 
be calculated in the incomplete Paschen- 
Back effect regime, which requires to di- 
agonalize the total Hamiltonian (see, e.g., 
Condon & Shortley 1970). Interestingly, 
as we shall see below, the proper calcula- 
tion of the Zeeman components results in 
significant differences in the Stokes profiles 
that emerge from a magnetized stellar at- 
mosphere. Furthermore, such differences 
do not vanish in the limit of weak fields. 

The outhne of this paper is as follows. 
Section 2 presents a brief review of the 
Paschen-Back effect theory, illustrating the 
range of magnetic strengths for which the 
J-levels of the 2'^P2,i,o term cross. Section 
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3 shows a detailed comparison of the re- 
sults that we have obtained assuming lin- 
ear Zeeman splitting (LZS) or incomplete 
Paschen-Back splittings (IPBS), including 
an analytical demonstration that explains 
our numerical radiative transfer results in 
a Milne-Eddington model. Finally, Section 

4 gives our concluding remarks. 

2. The Paschen-Back effect 

According to the theory of the Zee- 
man effect (see, e.g., Condon & Short- 
ley 1970) the corrections to the degener- 
ate energy Ej of any atomic level of to- 
tal angular momentum quantum number 
J are obtained by diagonalizing the ma- 
trix < JM\Hb\JM' >, where Hb is the 
magnetic Hamiltonian, given by: 

Hb^I^o{L + 2S)B, (1) 

with Ho = 9.27 x lO^^^ erg G^^ being the 
Bohr magneton. By choosing the quanti- 
zation z-axis of total angular momentum 
along the magnetic field vector, one finds: 

< JM\Hb\JM' >= iioBgM5MM' (2) 

where g is the Lande factor. This equation 
shows that, to first order of perturbation 
theory, any atomic level of total angular 
momentum quantum number J is split by 
the action of a magnetic field into (2J + 1) 
equally spaced sublevels, the splitting be- 
ing proportional to the Lande factor g and 
to the magnetic field strength. It is impor- 
tant to note that a state of the form | JM > 
is an eigenvector of the total Hamiltonian 
{Hq + Hb), with eigenvalue Ej + iioBgM, 
only when the z-axis of the reference sys- 
tem points in the direction of the magnetic 
field. 



The previous well-known results of first- 
order perturbation theory are correct only 
if the splitting produced by the magnetic 
field on a J— level is small compared to 
the energy separation between the differ- 
ent J— levels of the (L, S) term under con- 
sideration. In other words, the standard 
theory of the Zeeman effect is valid only in 
the limit of "weak" magnetic fields. Here, 
"weak" means that the coupling of either 
the spin or the orbital angular momen- 
tum to the magnetic field is weaker than 
the coupling between the spin and the or- 
bital angular momentum (the spin-orbit 
coupling) . This is the so-called Zeeman ef- 
fect regime. 

In the opposite limit, the magnetic field 
is so "strong" that the spin-orbit interac- 
tion can be considered perturbation 
compared to the magnetic interaction. In 
this case the magnetic field dissolves the 

— * 

fine structure coupling - that is, L and 

— * 

S are practically decoupled and precess 
independently around B. Therefore, the 
quantum number J loses its meaning. In 
this so-called complete Paschen-Back effect 
regime the magnetic Hamiltonian is diag- 
onal on the basis \LSMlMs >, and the 
eigenvalues are given by: 

< LSMlMs\Hb\LSMlMs >= iioB{Ml+2Ms) . 

(3) 

In this "strong" field regime, the term 
(L, S) splits into a number of components, 
each of which corresponds to particular 
values of {Ml + 2Ms). In the presence 
of a magnetic field, the corrections to the 
energy of the state \LSMlMs > due to 
the spin-orbit interaction has the form < 

— * — * 

AL ■ S >— AMlMs, where A is a constant 
with dimensions of energy. 

It is of interest to note that, since the 
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spin-orbit coupling increases rapidly with 
increasing nuclear charge, the conditions 
for a "strong" field are met at a much lower 
field with light atoms (like helium) than 
with heavy atoms. In fact, as shown in 
Fig 1, the levels with J — 2 and J — 1 
of the upper term of the He I 10830 A 
multiplet cross for magnetic strengths be- 
tween 400 G and 1500 G, approximately. 
This level-crossing regime corresponds to 
the incomplete Paschen-Back effect regime, 
in which the energy eigenvectors are grad- 
ually evolving from the form \LSJM > 
to the form \LSMlMs > as the magnetic 
field increases. This range between the 
limiting cases of "weak" fields (Zeeman 
effect regime) and "strong" fields (com- 
plete Paschen-Back regime) is more diffi- 
cult to analyze since it implies the eval- 
uation of matrix elements of the form < 
JM\Hb\J'M' >, which (in general) can 
only be done numerically. To this aim, 
we have developed a numerical code capa- 
ble of computing, for a magnetic field of 
arbitrary intensity, the strenghts and the 
splittings of the a and tt components in a 
line multiplet arising from the transition 
between two (L, S) terms. This code is 
very similar to a previous one that solves 
the problem of the intermediate Paschen- 
Back effect in a line with hyperfine struc- 
ture (Landi Degl'Innocenti 1978). 

3. Results 

We have developed a numerical code 
for the synthesis and inversion of Stokes 
profiles in a Milne-Eddington atmosphere, 
considering IPBS as described in §2. This 
code has been built with the aim of ana- 
lyzing spectropolarimetric observations of 
solar active regions, but we shall use it here 



to explore the effects of IPBS on the 10830 
polarization profiles. 

We consider a reference model atmo- 
sphere and calculate the multiplet com- 
ponent positions and strengths for various 
magnetic fields (ranging from to 3 kG). 
The Zeeman components are calculated ac- 
counting for IPBS and then compared to 
the usual approximation of LZS. The most 
relevant parameters of the chosen Milne- 
Eddington model are -q^ — A, AA^) = 
105 mA, Bq = 0.8, Bi = 0.2, a = 0.8 (sec, 
e.g., Landi Degl'Innocenti 1992 for defini- 
tions of the Milne-Eddington parameters). 
Sample profiles are depicted in Fig 2 for a 
500 G field with an inclination of 50° from 
the line of sight. There arc a few obvi- 
ous differences between the IPBS and LZS 
Stokes V profiles: a) The blue component 
of the multiplet, Trl, shows a net asymme- 
try with a global offset towards negative 
values;^ b) A small non-Doppler blueshift 
of the same component (which is actually 
related to the asymmetry); c) A reduced 
amplitude of the red and blue lobes of the 
red component of the multiplet (Tr2 and 
Tr3). The linear polarization profiles ex- 
hibit also some small asymmetry. These 
differences are due to IPBS effects on the 
Zeeman components of the transitions, dis- 
cussed in §3.1 below. 

In order to compare the IPB results 
with the LZS ones, one needs the rela- 
tive strengths and effective Landc factors 
of the three components. These are listed 
in Table 1. Note that some of the val- 
ues do not coincide with those used by 

^The sign of this offset depends on the polarity of 
the profile. It is produced by an imbalance in the 
weights of the a components of this transition, 
with the red component being stronger. 
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Table 1 

Atomic parameters used in the calculations 



A (A) 


Lower term 


Upper term 


9eff 


Relative strength 


10829.09 






2.00 


0.111 


10830.25 




^Pi 


1.75 


0.333 


10830.34 




=^P2 


1.25 


0.556 




Field (G) 

Fig. 1. — Energy diagram for the upper 
levels of Tr2 and Tr3 (the "^Pi and lev- 
els of the multiplet). Solid: IPBS. Dashed: 
LZS. Note the numerous level crossings be- 
tween 400 and 1500 G, approximately. En- 
ergies are in cm~^ with respect to the J=0 
level. 



other authors in previous works. For ex- 
ample, concerning the LZS calculations the 
effective Lande factor (geff) of Tr3 is 1.25, 
whereas Riicdi ct al. (1995) and Lagg et al. 
(2004) take it to be 0.875. The relative 
strengths^ (/i) of Trl, Tr2 and Tr3 are here 
0.111, 0.333 and 0.556, respectively. Our fi 
agree with those of Riiedi et al. (1995), 
but not with Lagg et al. (2004) (these 
authors quote relative oscillator strengths, 
but the conversion is straightforward) . It is 
not completely clear to us what the source 
of the discrepancy is. For the geff we 
used its standard definition (e.g., Landi 
Degl'Innocenti 1992) applied to the quan- 
tum numbers shown in the table. The 
fi are calculated according to the usual 
formulae (e.g., Radzig & Smirnov 1985). 
Thus, both of these quantities are derived 
straightforwardly from well-known atomic 
parameters. As we discuss below, it is not 
possible to reproduce the observations us- 
ing the LZS approximation. It is conceiv- 
able that other authors had to modify their 
atomic parameters in an attempt to com- 
pensate for IPBS effects. This would allow 

^Note that the symbol fi, which is used here and 
in other previous works on this multiplet, refers to 
the relative strengths of the components. These 
are not the "oscillator strengths" of the transi- 
tions, which are also commonly denoted by 
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them to fit observed profiles with those cal- 
culated using LZS. For example, the ratio 
of Stokes V amplitudes of the blue and red 
components calculated with LZS is signifi- 
cantly larger than it should be. Artificially 
reducing the g^ff of Tr3 reduces the ampli- 
tude of the red component without affect- 
ing the intensity profiles, thus bringing the 
profiles closer to the observations. 

3.1. Zeeman components 

The positions (with respect to the zero- 
field wavelength) and strengths of the var- 
ious TT and a components of the multiplet 
are shown in Fig 3. Tables 2 to 4 list the 
actual IPBS values at 500 G, which are also 
useful to match the curves in the left and 
right panels. Note that the positions of the 
Trl components arc not significantly dis- 
turbed by the presence of the other levels 
(J=l and J=2) in the field range depicted 
in the plots. However, the strength of the 
(T components is indeed perturbed, intro- 
ducing an asymmetry between the left- and 
right-handed polarizations. For our discus- 
sions in §3.2 below it is important to note 
that this asymmetry increases linearly with 
the magnetic field. 

The positions of the Tr2 and Tr3 com- 
ponents feel a much stronger perturbation 
for fields stronger than ~500 G. This is 
due to the close proximity in the energy 
diagram of the upper levels of these tran- 
sitions (the J —1 and J = 2 levels) . Sim- 
ilarly to the case of Trl, the component 
strengths of Tr2 and Tr3 become asym- 
metric. The asymmetry is rather linear for 
weak fields. As the field increases above 
~500 G, however, the behavior of the com- 
ponent strengths becomes non-linear. 

Finally, it is interesting to note the pres- 



ence of a TT component in Tr2 that does not 
exist in the LZS approximation. This new 
component appears for non-zero fields in 
IPBS and becomes more important as the 
field increases. 

The asymmetries between left- and 
right-handed a components in all three 
transitions, are important. They ulti- 
mately manifest themselves as asymme- 
tries in the profiles and lead to net circular 
polarization even in the absence of velocity 
or magnetic field gradients. 

3.2. Weak-field asymptotic behav- 
ior 

Fig 3 shows that the Zeeman compo- 
nents of the various transitions tend to- 
wards the LZS values as the magnetic 
strength goes to zero. One may thus ex- 
pect that the IPBS profiles (plotted in the 
top row) would become similar to the LZS 
profiles in the weak field limit. However, 
this is not the case. The profiles obtained 
from our calculations maintain their shape 
unchanged down to the lower limit of our 
simulation (0.1 G). This apparent paradox 
has a simple, albeit slightly counterintu- 
itive, explanation. The Stokes V profiles 
can be described conceptually as the sub- 
traction of two bell-type functions, simi- 
lar to Gaussians or Lorentzians. In the 
weak field limit the separation between 
these functions is much smaller than their 
width, and the resulting Stokes V profile is 
essentially a small residual of the subtrac- 
tion. A slight asymmetry between the two 
functions will show more conspicuously the 
closer together they are. Therefore, reduc- 
ing the field results in "amplifying" the 
red/blue asymmetries. Let us analyze this 
in some detail considering a very simple 
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Fig. 2. — Simulated Stokes profiles of the 10830 A multiplet, assuming a 500 G field inclined 
50° with respect to the line of sight. The profiles are normalized to the quiet sun continuum 
intensity (Ic). 



Table 2 

Component positions and strengths for TrI (500 G) 



lype of 


Position 


Strength 


component 


(Doppler units) 




a+ 


0.524 


0.106 


TT 


0.004 


0.111 


(7~ 


-0.516 


0.116 
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Tri . Component positions 




500 1000 1500 2000 
Field (G) 



Tr2. Component positions 




500 1000 1500 2000 

Field (G) 



Tr3. Component positions 




500 1000 1500 2000 
Field (G) 



0.20 
0.15 
0.10 

0.05 
0.00 



Tri. Component strengths 



500 1000 1500 2000 
Field (G) 



Tr2. Component strengths 




500 1000 1500 2000 
Field (G) 

Tr3. Component strengths 




500 1000 1500 2000 
Field (G) 



Fig. 3. — Positions and strenghts of the Zeeman components as a function of the magnetic 
field. In all panels the solid lines represent calculations considering IPBS while the dashed 
lines represent those using the LZS approximation. Note that, in the general IPBS case, the 
Zeeman components exhibit asymmetric displacements and strengths. Also notice that the 
Zeeman IPBS components have strengths that depend on the magnetic field. 



8 



Table 3 

Component positions and strengths for Tr2 (500 G) 



Type of 


Position 


Strength 


component 


(Doppler units) 




(7+ 


0.542 


0.136 


a+ 


0.410 


0.118 


TT 


0.149 


0.118 


TT 


0.022 


0.007 


TT 


-0.111 


0.215 




-0.371 


0.215 


a" 


-0.499 


0.190 



Table 4 

Component positions and strengths for Tr3 (500 G) 



Type of 


Position 


Strength 


component 


(Doppler units) 




(7+ 


0.495 


0.091 


(7+ 


0.371 


0.215 


(7+ 


0.260 


0.333 


TT 


0.111 


0.215 


TT 


-0.025 


0.215 


TT 


-0.149 


0.118 


a~ 


-0.260 


0.333 


a~ 


-0.410 


0.118 




-0.546 


0.027 
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scenario. Suppose that we have a Stokes V 
profile that is the result of subtracting two 
Gaussians separated by a small distance 
2Ax: 

V{x) = >le-^^-^^)' - Se-(^+^^)' . (4) 

The zero-crossing wavelength (xq) of V{x) 
is given by the condition: 

which, taking logarithms on both sides, 
leads to: 

log{B/A) 

If the asymmetry is small {B/A ~ 1), a 
Taylor expansion of the logarithm yields: 



Xo 



B-A 
ABAx 



(7) 



For weak fields, the asymmetry B — A de- 
pends linearly on the field (Ax). We then 
conclude that: 



Xq oc 



1 



i.e., the zero-crossing wavelength is shifted 
and does not depend on the magnetic field. 
In the derivation above we have been con- 
cerned only with strength asymmetries, ne- 
glecting the effect of asymmetric compo- 
nent positions. If we denote by Axa and 
Axb the (different) component positions, 
a similar straightforward reasoning can be 
conducted to conclude that: 

B-A 1,, , , 

where Ax = {Axa + Axb)/'^- It is evi- 
dent from Eq (9) that position asymme- 
tries have only a second-order effect on xq, 
since {Ax a — Axb) is proportional to Ax. 



We verified this property with our simula- 
tions. 

Consider now the area asymmetry: 



A 



(10) 



In the weak-field limit, the denominator of 
Eq (10) is directly proportional to the field. 
Thus: 



^-^1 V{x)dx, (11) 



where c is a constant. Expanding the ex- 
ponentials in Eq (4) in a Taylor series and 
truncating to first order, we obtain that for 
Ax < 1: 

^-{x-Axf ^ g-x2 ^ 2x0-""' Ax . (12) 
Thus: 

V{x) ~ (^ - B)e-^' + 2a;e^^'Aa;(yl + B) . 

(13) 

Inserting this into Eq (11) and taking into 

account the odd parity of the second term 
in Eq (13), it is evident that: 



A 



TT- 



A-B 
Ax 



(14) 



Therefore, the net area asymmetry A re- 
mains constant even for fields approach- 
ing zero. All of the asymptotic properties 
of Stokes V derived in this section (con- 
stancy of the zero-crossing wavelength, net 
asymmetry, relative importance of weight 
and position asymmetries) have been con- 
firmed by further modeling using our nu- 
merical code. It is important to note that 
the results we present here do not vio- 
late the principle of spectroscopic stability, 
since the Zeeman components (and their 
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center of gravity) approach asymptotically 
the LZS values. 

The asymptotic values that we obtained 
for Trl are a;o=340 m s^^ and |^|=0.11 
(with the red lobe having a larger area). It 
must be noted that this value of A is influ- 
enced by the blend with the blue wing of 
the red component. We estimated the ef- 
fect of this blend by calculating A for the 
LZS profile (which should be zero), obtain- 
ing |^LZ5|=0.04. In this case it is the blue 
lobe that has larger area. Therefore, the 
presence of the blend has the effect of re- 
ducing the of the IPBS profile. The 
actual value of xq is significantly model- 
dependent. In particular, it is very sensi- 
tive to the Doppler width of the line. Dou- 
bling this parameter in our model, resulted 
in a factor 4 increase of Xq. For the red 
component of the multiplet (Tr2 and Tr3), 
we obtained xo—95 m s~^ and |^|=0.05. 
Here, as opposed to Trl, the zero-crossing 
is redshifted and the blue lobe has a larger 
area. The dependence of Xq and |^| for 
both lines are plotted in Figs 4 and 5. 

4. Conclusions 

The calculations reported in this paper 
show that IPBS effects need to be taken 
into account for proper diagnostics based 
on the He I multiplet at 10830 A. This is 
a very interesting set of lines, with some 
unique properties that make it an ideal tool 
for the exploration of chromospheric and 
coronal magnetism. Care must be taken, 
however, to consider all the relevant phys- 
ical processes that intervene in the forma- 
tion of the lines. 

The (still few) works existing in the 
literature regarding Stokes observations 
of these fines do not seem to agree on 



500 r 




oL ^ ^ ^ ^ ^ ^ ^ ^^^^ 

1 10 100 1000 

Mognetic Field (G) 

Fig. 4. — Relative Stokes V zero-crossing 
shift between IPBS and LZS profiles (ab- 
sohitc value), function of the field 

strength. Stars: blue component (Trl). 
Diamonds: red component (Tr2-|-Tr3). 
The Trl shift is towards the blue, whereas 
the Tr2-|-Tr3 shift is towards the red. 




0.02 - 

0.00 [ ^ ^ ^ ^ ^ ^ ^ 

1 10 100 1000 

Mognetic Field (G) 

Fig. 5. — Area asymmetry (|^|) of the 
IPBS profiles as a function of the field 
strength. Stars: blue component (Trl). 
Diamonds: red component (Tr2-|-Tr3). 
Trl has a larger red lobe, whereas Tr2-|-Tr3 
has a larger blue lobe. 
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such fundamental parameters as: the rest 
wavelength of the Trl line; the effective 
Lande factor of the Tr3 line; the relative 
strengths of the transitions (compare the 
values given in Riiedi et al. 1995; Lagg 
et al. 2004; this work). These are basic pa- 
rameters that should be well known, either 
from laboratory experiments or from very 
simple calculations. We believe that, at 
least partly, the existing confusion might 
have its origin in neglecting IPBS effects. 
We have shown that IPBS is responsible for 
a non-Doppler blueshift of Trl, systematic 
profile asymmetries and a reduced ampli- 
tude of Stokes V in the Tr2+Tr3 blend. 
As an illustrative example, consider the 
simple line-ratio experiment represented 
in Fig 6. For a fixed model atmosphere, 
the ratio of Stokes V amplitudes in the red 
(Tr2-|-Tr3) and the blue (Trl) components 
of the multiplet depends smoothly on the 
magnetic field. It is clear from the figure 
that neglecting IPBS leads to an erroneous 
estimate of the field. Interestingly, the 
IPBS curve is steeper than the LZS one, 
indicating that it is more sensitive to the 
field. 

Some IPBS signatures described in this 
paper are probably present in published 
observations. If one compares the Trl 
Stokes V profile depicted in our Fig 2 
(upper right panel) with those of Lagg 
et al. (2004) (Figs 4 and 5), the resem- 
blance is striking. More specifically, the 
IPBS blueshift and the net offset towards 
negative values both become immediately 
apparent. Profile asymmetries are also 
present in their observations. However, it 
is not possible to assess whether these arise 
from IPBS or from velocity or magnetic 
field gradients. 
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Fig. 6. — Ratio of the Stokes V amplitudes 
of the red (Tr2+Tr3) and blue (Trl) com- 
ponents of the multiplet, as a function of 
the magnetic field for a given model. Solid: 
IPBS. Dashed: LZS. 
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